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Multiplication of a inatrix by its transpose^ either 

m rn 

AA =C or A A= C* 

Entire program 

If ^ <6k max i = 160 

If j < 128 max i = So 

10 decimal places 

The routine will read a matrix. A.., and form either 

AA = C or A A = C* as specified by the parameter tape. 

Featiires of the routine are; 

1. The input matrix A, can be rescaled by any power of 
ten prior to multiplication, (provided the elements re- 
main within the range of the machine -1 < a^ < +l) 
enabling the user to adjust the scaling on C either 

to prevent overflow or for other reasons. The machine 
also checks for overflow of elements of C- (See section 
on stops and error diagnosis). 

2. The output matrix, C or C*, will be symmetric. The 
parameter tape specifies four output options? 
triangular in row or column form, 

or square in row or column form. 
5o If overflow should occur, proper switching allows a new 
scale parameter to be inserted and the problem redone 

without need to re-read data. 

Stops 

1. ^fester tape 5^08? 

2 . Parameter tape 240K6 

3. Matrix A tape 2^087 

At stop 2^087, a new parameter tape can be read by raising 
black switch. To change the scaler only, move white switch 
up and down, then raise black switch. (See section on stops 
and error diagnosis). 
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PARAMETER TAPE PREPARATION: 

There are f oiir parameters with terminating symbols and 

they must appear in the following order: 

dJ tF xL +sS. 

dJ: d is the number of decimal places desired in the 

output matrix. 

T 
tF: When t = 0, A A will be formed 

T 
when t = 1, AA will be formed 

xL: X determines the form of the output matrix C or C*« 
If X =0, the output will be in triangular form by 
columns punched with a carriage return after each 
element and two carriage returns at the end of each 
row. 

If X =1, the output will be in triangular form 
punched by rows across the page| two carriage returns 
separate rows, and if the row is too wide for the 
page, appropriate carriage returns are inserted. 
If X =2, the output will be square by columns with 
a carriage return after each element. An N and 
two carriage returns are punched after each row. 
If X = 3, the output will be in square form by rows 
across the page with an N and two carriage returns 
after each row; if the row is too wide for the page, 
appropriate carriage returns are inserted. 

+sS: +s or -s is the integer exponent of ten, the 
or 

oS: scaling constant. Each element' of A is multiplied 

+s 
by 10-- prior to multiplication. If no rescaling 

on A is desired, this should appear as oS. (NO 

SIGN should be indicated). 
DATA TAPE PREPARATION: 

The elements of matrix A are pimched by rows as signed 
fractions with an N at the end of each row. If an F is 
punched instead of an N, the computer will stop on either 
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500K8 or 300SS. By raising black svitch reading of A 
is resumed. After the final N, a J must be punched to 
indicate the end of the matrix. 
STOP AND ERROR DIAGNOSIS: 

DIAGNOSIS 
Master tape has read correctly; black 



STOPS 
5^87 



LOG. 



2kOK6 
5OOK8 
3OOSS 
24087 



0Ki4- 
0K8 
OSK 
13J 



switch to read the parameter tape. 

End of parameter tape; black switch to read 

in data; white switch to skip data input. 

F termination, first row; black switch to 

continue . 

F termination J subsequent rows; black switch 

to continue. 

End of problem 

(a) Black switch to begin a new problem 
with a different parameter tape. 

(b) To rescale only raise white switch up 
and down and then raise black switch to 
read in new scaling parameter. 

(c) If just new data is to be read in using 
former parameters, move white switch up 
and down twice, then black switch to 
read in new data. 

Number of elements in rows are not equal. 
Check data tape. 

FF002 11F,171 Overflow in an element of C. To read in a 

003 122-175 

new scaling parameter -sS only, move white 

switch up and down; then to skip re-reading 

the data white switch again. 

m 

Sum check failure in matrix A or A . White 
switch three times; re-insert the data, and 
then black switch to start program. 
Master tape sum check failure, re-read 
master tape. 



FFOOl 



OSF 



FFOOi^ 
005 



15K 
13J 



FFOO6 5F5 
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T 



DURATION IN SECONDS: AA"" square form 



i (.02d + .002 j + .12) + i (.005dJ + .Ol8j + .05) +40 
triangular form 

i (.Old + .002 j + .07) + i (.005dj + .Ol8j + .05) +4o 
A A square form 

J (.02d + .0021 + .12) + J (.005di + .0181 + .015) +^ 

triangular form 

j (.Old + .002i + .07) + J(.005di +..0l8i + .015) +^0 



nj 



DAT E June 27, I9S0 
PROGRAMMED BY 



APPROVED BY 






LOCATION 




(BDER 


NOTES 
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KSL 5.02 


Abs. 


Rel. 


Sym. ■" 


003K 














5 





' 


OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 


d 

T T 
0, 1 A-^A, AA"" 

0, 1, 2, 3 form of 

J 

65 or 129 

Sum of elements of 

I 


output 
A 




10 


10 




OOF OOF 
85IIF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
OOF OOF 
85IIF OOF 
OOF 00(1) 
OOF OOF 


2560 - 129 or 65 

i 

J 

k 

Sum of elements of 

Sum of elements of 

[A^,Ior[A^.] 

simi box 


A 
B 




20 






OOF OOF 


tally of row elements 










OOF OOF 


65/d + 2 or 1 












OOF OOF 


Scaling factor 








20 




OOF OOF 
7J22F 40600F 
S5F i^0600F 


21(N12) 












OOK 








26 





(1) 
(2) 
(6) 
(9) 


OOF OOIF 
OOF 002F 
OOF oo6f 
OOF OO9F 








30 




(10) 

(65) 

(129) 

(600) 

(2560) 


OOF OOIOF 
OOF OO65F 
OOF OOI29F 
OOF OO6OOF 
OOF OO256OF 











LOCATION 


ORDER 
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KBL 5.02 




Abs. 


Rel. 


Sym. 


80F OOF 








35 


9 


(-1) 








(Dl) 


OOIF OOIF 












(TOD) 


86nF OOF 












(FRD) 


8511 F OOF 












(end) 


i4-0(TRC) L56OOF 
OOK 








ko 




(N12) 


OOK 








79 




(P16) 


OOK 








135 





(IN) 
(PAR) 


I92F k02iF 
kll' QlkY 
LO(lO) 36(STO) 
Li+(10) 5OF 
7^^(10) S5F 


Read first three parameters 






llfO 




(STO) 
(STA) 


kOF 22(PAR) 
if2(STA) 92575F 
L5F kOF 
L52F L^2F 
k02F 36(PAR) 
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(sc) 

(INK) 
(IND) 


I^-IF 8lifF 

LO(io) 36(im:) 

ifl22F 22(DATA) 
1+OF ^IIF 
8lifF LO(lO) 


Read in scaling parameter 






150 


20 


(IKE) 
(INH) 

(inf) 


32(IKE) Lif(lO) 
50IF 7i4-(lO) 
S5F if OIF 
26(IKD) T.11F 

k-02F 50(10) 
F52F i4-02F 
32(IWF) 75(10) 
26(INH) S5F 
if02F L3F 


Compute 10 





LOCATION 


ORDER 
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.02 


T 


Abs. 


Rel. 


Sym. 


36(ING) kW 








159 


25 






160 


26 


(ING) 


50FL5(l) 
662F S5F 
40??F 22(DATA) 
L52F lA(-l) 


Compute 10~^ 








30 


(data) 

(DAT) 

(tes) 


i|-022F 2^(DAT) 
22(PA) 26(DAT) 

50600F 50(dat) 

26(N12) L0(2) 
30(TES) i4-l6F 
L521(N12) 4023F 


Black switch to read in datt 
White switch to skip data 
Read first row- 
Stop for F at end of first row 






170 






1020F L0(600) 
i|-26F L56F 
L0(65) 32(DRU) 


Store J 








ko 


(dru) 

(SCH) 


L5(65) ^07F 
22l(DRU) L5(129) 
i^OTF L5(6) 
i^2(SCKB) L5(l) 
i<-09F ifl8F 
50600F 50(SCH) 
26(SCK) L5F 


Store 65 or 129 into L Tj] 
Set I = 1 






180 


50 


(ina) 
(tis) 


hoQF L5(2560) 
L07F 4010F 
L5(600) Lk6F 
i<-2(END) 26(TRA) 
506OOF 50(IKA) 
26(N12) L0(2) 
301(TIS) Lif(l) 
32(PA) L521(N12) 
L023F 40F 
L3F 32(INB) 


Store sum first row 
Store 2560 - 129 or 65 

Read in rest of data 
Stop on F 
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(INB) 
(INC) 


FFIF 26(lNC) 
506OOF 50(INC) 


FFl - wrong number of elements in row 
Sum row and add to sum check box 












26(SCK) L5F 











LOCATION 


ORDER 
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Abs, 


Rel. 


Sym. 


LkOF k08F 








193 


59 








60 


(TRA) 

(trb) 


F59F i<-09F 
L5(600) i|-2(TRB) 
L510F L47F 
il-OlOF lA(TOD) 
i]-0(TRC) L56OOF 


1 = 1 + 1 










(•mc) 


86IIF OOF 


Store rov starting at L 6OO, 


to drum 




200 


70 


■ 
(PA) 


F5(TRC) 40(TRC) 
F5(TRB) 1^2(TRB) 

lo(end) 36(ina) 

22(TRB) L55P 
0039F 36(pf) 
L53F 14(2) 
OO7F 40F 
L5(65) 66f 
S5F IO32F 
4021F 26i(PF) 


Form 65/d + 2 






210 


80 


(PF) 
(PG) 


l^[l) kOPlY 
L53F 0020F 
46(PRB) L3?1F 
4020F 92139F 
L54F LO(l) 
36(AT) 26(TA) 

OOK 


Store d into print routine 

Switch to Aa'^ or a'^A 
Drum transfer subrbutine 






216 





(DR) 
(DRD) 


K5F 42(DRA) 
46(DRB) kOF 
F5F 42(DRD) 
42(RET) L5F 


■ 






220 




(DRA) 

(drb) 

(DRF) 

(dre) 


46(DRF) L5F 
42{DRC) 26(DRB) 

LkY 42 ■'fin) ■ 
L5F 4o(DRE) 
85IIF OOF 


1 





LOCATION 



Abs, 



225 



Rel. 



231 



9 
10 



2if0 



250 



sym. 



(drc) 
(drk) 

(red) 
(fin) 

(SCK) 
(SCKA) 



(SCKB) 

(res) 
(sum) 

(at) 



ORDER 



10 



(ATJ) 



(ATA) 



(ATB) 



JOF kOF 
tojDREj L5/rEC) 

i^2Cdrc; loCfinJ 

36(DRE) 22F 
5OF 40F 

OOK 

K5F i^2(RES) 
1020F ^2(SCKA) 
i^2(SUM) i+lF 
L5F L4f 
1+OF F5(SCKA) 
1|2(SCKA) L5(SCKA) 
L0( :Ul.l) LOoF 
36(SCKA) 22F 
75F l4f 

OOK 

L5(6oo) l46f 

k22k-F if225F 
L518F ^2(dRK) 
L5(6) k2(SCKB) 

L5(frii) l4(256o) 

i^OllF L5(l) 
4012F ^^15F 
L5(l) ^OlifF 
L5(FRD) l4(2560) 
^017F L5(600) 
if2(ATP) 42l(ATP) 
i+2(ATQ) i<-2l(ATQ) 
506F 50(ATA) 
26 (DR) OO6OOF 
OOllF 4116F 
50600F 50(ATB) 
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CD Planted from (TA) or (AT) determines 
whether subroutine reads row or column 



Sum check of matrix rows subroutine 



Either 6 or 9 planted from (AT) or (TA) 



AA 



T 



Transfer row A to position starting L ^OOJ 



Sum elements in row A . ., 

" il 



LOCATION 




ORDER 
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Abs. 


Rel. 


Sym. 


26(SCK) L5F 




256 


16 










Li^l5F i|-015F 










L322F 32(ATS) 








(ATU) 


L522F 32(ATQ) 


Multiply rov "by scaling factor 


260 


20 


(ATP) 
(ATQ) 


1039F 7560OF 
S5F i+0600F 
F5(ATP) i4-2(ATP) 
F5l(ATP) i4-2l(ATP) 
L025F 32(ATS) 
26(ATU) 506OOF 
7J22F 14-0600F 
F5(ATQ) i<-2(ATQ) 
F5l(ATQ) 421(ATQ) 
L02ifF 32(A3S) 




270 


30 


(ATS) 


22(ATQ) L5(p) 
^(ATE) i^2(ATE) 








(ATC) 


506F 50(ATC) 
26(DR) OO8OOF 
0017F i<-119F 


Transfer rov A to position starting lJ6 






(atdO 


508OOF 50(ATD) 
26(SCK) L5F 
Li+l6F kOl6F 
L5(l) if013F 


Sum elements in rov A 






(ate) 


506OOF 7J8OOF 


Italtlply A A 


280 


ko 


(Atz) 

(ATI) 


kO¥ 36i(aty) 

L519F 32(ATZ) 

262(atz) l4f 

i^019F 26(ATW) 
Lk¥ 36(ATY) 
i^019F 26(ATW) 
FF2F 26(SC) 
L519F 32(ATX) 
LkY k019Y 


Start of test for overflov 






(ATX) 


26(atw) iAf 





LOCATION J 


ORDER 
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KSL 5 •02 


Abs. 


Rel. 


Sym. 


322 (ATX) i''J?'3F 






290 


50 










26 (SC) 1+019F 










(ATW) 


L513F L06f 
56 (PRA) L5(ATE) 
lA(Dl) i^O(ATE) 
F515F hQl3Y 


Print C^jj. 








(ATO) 


26 (ATE) L55F 
L0(2) 36(A) 
L51^F L012F 
361 (ATH) 221(A) 






500 


60 


(A) 


L5li4-F LO9F 
52 (ATH) I^ITF 
I^TF i+OlTF 
F51^F iK)l^F 










(ATH) 


22 (ATS) 927T0F 
92135F 92519F 
L32IF 4020F 
L312F LO9F 
32 (ATK) L511F 
IA7F 4011F 


Print N at end of row 
Reset row tally 




310 


TO 




F512F i4-012F 










(ATK) 


26(ATJ) L58F 
L015F i^OF 
L3F 32 (ATL) 


Start of sum checks 








(ATL) 


ff4f L58F 

LOI6F i+OF 
L3F 32 (ATM) 










(ATM) 


FF5F 2^4- (in) 
2U(SC) 2J+(DAT) 


End of problem - stop 2^4-087 








(PRA) 


5OF 50 (PRA) 






320 


80 


(P) 


26 (PR) 22 (ATO ) 
OO6OOF OO8OOF 







LOCATION: 


ORDER 
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khsi 


Rel. 


Sym, 


OOK 


T 
A-^A 








522 





(TA) 
(TAJ) 


L5(6oo) l49F 
4224F 4225F 
L5(6) l4(i) 
42(DRK) L5(9) 
42(SCKB) L5(FRD) 
l4(2560) 4011F 
L5(l) 4014f 
4115F L5(l) 




550 


10 




4015F L5(FRD) 
l4(2560) 4017F 
L5(600) 42 (tap) 
421 (tap) 42 (taq) 
421(TAQ) 42F 








(taa) 


509F 50 (taa) 
26 (dr) 00600F 
OOllF 4116f 


Transfer coliamn A., to position starting 






(tab) 


50600F 50 (tab) 


Sim elements in colxMn A.., 








26(SCK) L5F 




540 






l415F 4015F 
L522F 52 (tas) 


Multiply column by scaling factor 




20 


(tau) 


L522F 52 (TAQ) 








(tap) 


1059F 7560OF 
S5F 40600F 
F5(TAP) 42 (tap) 
F51(TAP) 421 (tap) 
L025F 52 (tas) 








(taq) 


26 (tau) 50600F 

7J22F 40600F 




550 


50 




F5(taq) 42 (taq) 

F51(TAQ) 421 (taq) 
L024F 52 (tas) 








(tas) 


22 (taq) L5(P) 





location 


ORDER 
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.02 


Abs ♦ 


Rel. 


Sym. 


1^6 (tae) 42 (tae) 






55^ 


32 










(TAG) 


509F 50 (tag) 
26 (dr) OO8OOF 

0017F 1+119F 


Transfer colimn A,, to position stai 


rtin 






(tad) 


50800F 50 (tad) 

26(SGK) L5F 


Sum elements of columns A. . 




360 






1A16F 4oi6f 

L5(l) h012F 








ko 


(tae) 

(TAZ) 

(tay) 


5060OF 7J8GOF 
il-OF 36i(TAY) 
L519F 32 (TAZ) 
262 (taz) iM 

iK)19F 26 (taw) 

iAf 36 (tay) 

lK)19F 26 (taw) 
FF2F 26 (so) 


Multiply A^j^ A^ . 

Start of test for overflow 




370 






L519F 32(TAX) 








50 


(tax) 
(taw) 

(tao) 


26(taw) iAf : 

322 (tax) FF3F 
26 (SC) 4019F 
L512F L09F 
36(prd) L5(TAE) 
iA(di) 14-0 (tae) 

F512F i4-012F 
26 (tae) L55F 


Print C* . 




580 






lo(2) 36(B) : 

L513F LOli^F 








60 


(B) 


361 (TAH) 221(b) 
L513F L06f 
32 (TAH ) L517F 






385 


63 




Li4-(l) il-OlTF 
F513F 4013F 






-, 




(tab) 


22(TAS) 9277OF 


Print N at end of row 





LOCATION 



Abs. 



588 



590 



Rel. 



66 



TO 



597 



iK)0 



996 



Sym. 



(tak) 

(PRD) 



(PR) 
(PC) 

(prb) 



(PRC) 
(PZ) 



ORDER 



92135F 92519F 
L521F k020F 
L5li4-F L06f 
52 (tak) L511F 
lA(l) it-OllF 
F5lil-F ifOli^F 
22 (TAJ ) 22(ATK) 
5OF 50 (PRD) 
26 (PR) 22(TA0) 

OOK 

K5F k2{FL) 
L520F 56 (PRC ) 
92963F L519F 

50F 50 (prb) 
26(Pl6) F52OF 
i<-020F 22 (PZ) 
9215IF L321F 
4020F 92519F 
26 (PC) 22F 

OO996K 
L3F 5^ (IN) 

ff6f 26 (in) 

K50T58F l^-90405F 
26L 26IN 
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Reset tally for rows 



Transfers to sum check routine in (AT) 



Print subroutine 



